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ABSTRACT 

Among active galactic nuclei, BL Lac objects show extreme properties that have been interpreted as the effect of relativistic beaming 
on the emission from a plasma jet oriented close to the line of sight. The Doppler amplification of the jet emission makes them ideal 
targets for studying jet physics. In particular, low-power BL Lacs (LPBL) are very interesting because they probe the jet formation 
and emission processes at the lowest levels of accretion. However, they are difficult to identify since their emission is swamped by the 
radiation from the host galaxy in most observing bands. 

In this paper we propose a new LPBL selection method based on the mid-infrared emission, in addition to the traditional optical 
indices. We considered the radio-selected sample of Best & Heckman (2012, MNRAS, 421, 1569) and cross-matched it with the 
WISE all-sky survey. In a new diagnostic plane including the W2-W3 color and the Dn(4000) index, LPBL are located in a region 
scarcely populated by other sources. By filtering objects with small emission line equivalent width, we isolated 36 LPBL candidates 
up to redshift 0.15. Their radio luminosity at 1.4 GHz spans the range logL^ = 39.2-41.5 [erg s“*]. 

Considering the completeness of our sample, we analyzed the BL Lac luminosity function (RLF), finding a dramatic paucity of LPBL 
with respect to the extrapolation of the RLF toward low power. This requires a break in the RLF located at log L, ~ 40.6 [erg s“*]. The 
consequent peak in the BL Lacs number density is possibly the manifestation of a minimum power required to launch a relativistic 
jet. 
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1. Introduction 


BL Lac objects (BL Lacs), together with flat-spectrum radio 
quasars (FSRQs), form a class of active galactic nuclei (AGNs) 
known as “blazars”. Like the other radio-loud AGNs, they are 
thought to be powered by a central supermassive black hole, 
surrounded by an accretion disk. Two plasma jets are ejected in 
the direction orthogonal to the disk. Among radio-loud AGNs, 
blazars show extreme properties, such as strong flux variability 
on the whole electromagne tic spectrum even on intra-day time 
scales ( A^gneret^ 1996I). variabl e optical and radio polariza¬ 
tion (lAller et al.l 1996l:~ SmithllT99^. and apparent sup erluminal 
motion of radio components dKellermann et al.ll2004ll . The ex¬ 
planation of this observing evidence requires that one of the rel¬ 
ativistic plasma jets be closely aligned with the line of sight, 
causin g Doppler beaming of the jet emission (lUrrv & Padovanil 
Il995l) . This makes blazars the ideal candidates for studying the 
physics of AGN jets. In particular, the least luminous blazars are 
extremely interesting because they can probe the jet formation 
and emission processes at the lowest levels of accretion. 


The classical distinction between BL Lacs and FSRQs is 
based on the strength of the emission lines in the optical spec¬ 
trum, with BL Lacs being identified by rest-frame equivalent 
widths (EW) less than 5 A (IStickel et al.lll99lh . The blazar spec¬ 
tral energy distribution (SED) is dominated by the non-thermal 
emission from the jet. It shows two bumps; the low-energy one 
(from radio to UV-X-rays) is due to synchrotron radiation, while 
the high-energy bump (X- and y-rays) is probably due to an 
inverse-Compton process off the s ame relativist ic electrons that 
produce the synchrotron photons (lKonigllll98lh . Depending on 
the frequency of the synchrotron peak Vp in the SED, BL Lacs 
have been divided into low-, intermediate-, and high-energy syn¬ 
chrotron peaked BL Lacs. The objects of the first class (named 


LBL or LSP) have Vp < 10'"^Hz. For those of the second class 
(IBL or ISP) 10'"^ < Vp < lO'^Hz, and for objects of the 
third class (HBL or HSP), Vp > lO'^Hz (lAbdo et aI.ll20Tol) . It 
has been claimed that blazar SEDs follo w a “blazar sequence” 
(iFossati et al.l 119981: iDonato et al.l 1200 ll) . linking the SEDs of 
FSRQs at the highest synchrotron peak luminosities Lp and low¬ 
est Vp to the SEDs of LBL, IBL, and Anally HBL, when pro¬ 
gressively decreasing Lp while increasing Vn. Thi s view has been 
criticized, for example, bv iGiommi et al.l(l2012ah . who underline 
the selection effects of flux-limited surveys and lack of redshift 
determination for the likely high-Lp-high-Vp sources. 

The difficulty of finding low-luminosity BL Lacs comes 
from the dilution of the jet component by the host-galaxy emis¬ 
sion in the optical and near-IR bands, as well as from the flux 
limit of the available surveys. As a result, their radio luminosity 
funct ion (RLF) is only determined above log L r ~ 40.5 [ergs '] 
(e.g.. IWo lter et al. 199^ t P adovani et al.ll2007h . The analysis of 
lMarcha& Caccianigal ( 20131) reaches lower luminosities when 
also including “Type 0” and weak-lined radio galaxies. Our aim 
is to And a substantial number of LPBL and to build their RLF 
down to a poorly studied regime of radio power. 

This paper is organized as follows. In Sect. 2 we present a 
new selection method for identifying BL Lacs that involves both 
optical and MIR properties. This leads to singling out 36 BL Lac 
candidates with redshift z < 0.15 and radio luminosities in the 
range logLi ~ 39 to 42 [ ergs“'] in Sect. 3. The various sources 
of incompleteness in our sample are considered in Sects. 4 and 5, 
while the resulting BL Lac RLF is derived and also discussed in 
the framework of the AGN unification scheme in Sect. 6. After 
a summary, our conclusions are drawn in Sect. 7. 

Throughout the paper we adopt a cosmology with Hq = 
70kms *Mpc"', Dm = 0.29, and Qa = 0.71. When we speak 
of radio luminosities Li, we mean vj^, where the monochromatic 
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Fig. 1. Left. Relation between the Dn(4000) index and the W2-W3 color measured from the sample of 21 065 bright nearby galaxies 
(see text for the selection criteria). For each bin in W2-W3, the large black dots denote the mean value and dispersion of the best 
Gaussian fit to the Dn(4000) distributions. Each spectroscopic class is highlighted with a color: red for lineless galaxies, black for 
LINERs, green for Seyferts, and blue for star-forming galaxies. The solid line represents the result of a toy model where a fixed 
fraction of a blue light component of increasing strength is reprocessed into MIR radiation. Right. Simulated Dn(4000) versus 
W2-W3 tracks obtained by adding a power-law component of increasing strength with a spectral index increasing from 0 (green) 
to 1 (cyan) with a step of 0.25 to the SED of a quiescent galaxy. Within each track we label the ratios of 0, 0.1, 0.3, 1, 3, and 10 
between the jet and galactic components. The black lines trace the iso-density contours of the distribution shown in the left panel; 
the red lines show the iso-densities contours for the objects with all lines rest frame EW < 5 A. 


luminosity is measured at Vi = 1.4 GHz. Moreover, we distin¬ 
guish between radio luminosities in cgs units, L^, and those in 
MKS units, R, . 


2. The selection method 


As mentioned, the peculiarity of BL Lacs is their non-thermal 
emission from a relativistic jet seen at a small angle with re¬ 
spect to the line of sight. As a consequence, their optical spec¬ 
trum is either featureless or only presents emission lines of small 
EW. This is because lines (when present) are isotropic, while the 
jet emission is highly beamed. Similarly, the jet dilutes the stel¬ 
lar absorption features. BL Lacs have been historically selected 
by se tting a limit of 5 A to the rest-frame EW of any emission 
line (IStickel et al.lll991l) . Alternatively, they have been recog¬ 
nized as objects with low Dn(4000) in detfl The Dn(4000) limit 
was set to 1.33 b ylStocke et al. ( 1991h. th en later relaxed to 1.67 
bv iMarcha et al.l (|1996l) and iPlotkin et al.l (l2010h FI Nevertheless, 
passive elliptical galaxies have Dn(4000)~2 (see below), so that 
we can expect that the Dn( 4000) constraint can be furt her re¬ 
laxed to find LPBL (see also lMarcha & Caccianigair2013l) . 


We propose to loosen the Dn(4000) limit but to include a 
further selection tool based on the mid-infrared emission (MIR). 
In fact, in active galaxies the MIR is also an isotropic quantity, 
owing to the reprocessing of (mainly) UV photons from the in¬ 
terstellar dust, while in BL Lacs it also receives a non-thermal 
contribution from the jet that alters the overall SED shape. As 


' Defined according to lBalogh et all lll999h as the ratio between the 
flux density measured on the “red” side of the Ca II break (4000^100 
A) and that on the “blue” side (3850-3950 A). 

^ Various authors use diiferent definitions of the Ca II break strength. 
Our parameterization is related to thos e found in the literature as 
Dn(4000) = 1/(1 - C) dLandt et al.ll2002ll . 


we discuss in detail below, the resulting MIR properties enable 
us to separate BL Lacs from other classes of active galaxies. 

The WISE satellite provides us with an all-sky survey in the 
four bands Wl, W2, W3, an d W4, which ^e cen tered at 3.4,4.6, 
12, and 22 /rm, respectively (IWright et al.l20ldl) . The best-suited 
color for our purposes is obtained from the W2 and W3 bands, 
because it has high wavelength leverage and level of complete¬ 
ness. 

Before we can proceed to use the MIR observations for the 
BL Lacs selection, we must explore the relationship between the 
optical and MIR properties of galaxies in general. Eor this pur¬ 
pose, we consider the 927 552 SDSS DR 7 spectra of the 818 333 
galaxies (MPA-JHU sample hereafter) in the value-added spec¬ 
troscopic catalogs produced by the group from the Max Planck 
Institute for Astrophysics, and The Johns Hopkins University, 
and available at http://www.mpa-qarching .mpg. de/SDSS/ 
dBrinchmann et al.ll2004t iTremonti et alJl2004 ). 

Our aim is to detect the faintest BL Lacs, so we started to 
analyze galaxies up to z = 0.1 and considered only massive ob¬ 
jects, since they are the typical BL Lac hosts. We then filtered 
objects with high-quality data with the following criteria: 

i) redshift z < 0.1, 

ii) error on Dn(4000) < 0.02, 

iii) signal-to-noise (S/N) in both W2 and W3 bands > 10, 

iv) stellar velocity dispersion > 100 kms *. 

In Pig.[T](left panel) we compare the W2-W3 color with the 
Dn(4000) index for the 21 065 selected galaxies. The presence 
of a clear sequence emerges, with the Dn(4000) value decreasing 
with increasing W2-W3. In each W2-W3 bin, we estimated 
the mean value and dispersion of the best Gaussian fit to the 
Dn(4000) distribution. 

Taking advantage of the SDSS spectra, we can consider the 
various classes of galaxies separately. They are located in dif¬ 
ferent areas of this diagram: the passive galaxies are clustered at 
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Fig. 2. Left. Relation between the Dn(4000) and H(5 indices from the SDSS subsample of 21 065 galaxies. For each bin in Dn(4000), 
the black dots indicate the mean value and dispersion of the best Gaussian fit to the Hb distribution. Right. Hb versus Dn(4000) 
for the outliers of the Dn(4000) versus W2-W3 relation (black dots). The red dots label, for each bin in Dn(4000), the mean value 
and dispersion of the best Gaussian fit to the Hd distribution for the overall SDSS galaxies population. The straight blue line is the 
expected locus of BL Lacs, obtained by varying the ratio between jet and galactic emission (see text for details). 


the top left, with Dn(4000) ~ 2 and W2-W3 ~ 0.5. Proceeding 
toward the opposite corner of the diagram we find LINERs, most 
of them having Dn(4000) ~ L5-2.0 and 1 < W2-W3 < 3, and 
then Seyfert galaxies with Dn(4000) ~ 1.5 and 3 < W2-W3 < 4. 
In the bottom right hand corner, we find star-forming galaxies, 
with low Dn(4000) values (~ 1.2-1.3) and large excesses in the 
W3 band. Therefore, all spectroscopic classes of galaxies follow 
the same overall sequence possibly because, regardless of the 
origin of the blue/UV light, the same mechanism of reprocess¬ 
ing into the MIR is operating. 

To support this interpretation further we built a very sim¬ 
ple model. We added a growing contribution of blue light to the 
SED of a passive galaxy and adopted a constant reprocessing ef¬ 
ficiency into MIR emission. The track obtained is shown in Eig. 
[1] (left panel) and follows the central values of the distributions 
remarkably well. 

It is important to note that the same sequence is followed by 
the objects with low values of the emission lines EW, the sources 
among which we expect, in principle, to find BL Lacs. By set¬ 
ting a limit to the lines rest frame of EW < 5 A, we obtain the 
red contours overplotted in Eig. [1] The distribution of this sub¬ 
sample follows that of the general galaxies population closely 
up to W2-W3 ~ 3. The low EW objects only avoid the region 
of largest mid-infrared colors, which are typical of star-forming 
galaxies. 

We can now explore the expected behavior of BL Lac ob¬ 
jects in the same diagram. This was simulated by adding to the 
SED of a quiescent galaxy a jet component Tjet. Eor the quies¬ 
cent galaxy we use a Dn(4000)=L98 and W2-W3=0.5, i.e. the 
coordinates of the upper left average point in Eig. [1] while for 
the jet component we adopt a power-law spectral shape in the 
MIR-optical frequency range of the form Tjet = Fq When 
Fjet increases, the W2-W3 color grows, owing to the larger rela¬ 
tive contribution of Tjet in the W3 band than in the W2 band with 
respect to the galactic emission. At the same time, an excess of 
blue light emerges, thus reducing the Dn(4000) index (see Eig. 
[T] right panel). Depending on the value of a, the expected lo¬ 
cus of BL Lacs varies. The curves obtained with different values 
of a can be overplotted onto the relation between W2-W3 and 


Dn(4000) obtained above. All the simulated tracks are steeper 
than the phenomenological trend traced by the sample consid¬ 
ered above. As a result, the expected BL Lac tracks exit the re¬ 
gion where the SDSS galaxies are located. The transition into 
the “forbidden zone” occurs at Dn(4000) ~ 1.8 for the flatter 
BL Lac SED and at lower values for softer SEDs. The visual 
inspe ction of broad-band SEDs of various cla sses of BL Lacs 
(e.g.. lGiommi et ani2012bHRaiteri et al.ll2014l) indicates that an 
optical-to-MIR index a ~ 0.5 is typical of HBL, while a ~ I 
can be associated to LBL. 

In other words, the presence of a BL Lac nucleus produces 
a substantial excess of blue light that corresponds to a smaller 
MIR excess than in the case of active or star-forming galaxies. 
We conclude that this diagnostic plane can be used to isolate BL 
Lac objects. 

2.1. Outliers in the Dn(4000) versus \N2-\N3 plane 

Although SDSS galaxies define a clear locus in the Dn(4000) 
versus W2-W3 plane, several outliers are present in the region 
where we expect to find BL Lacs. These objects thus deserve 
further investigation in order to ensure the purity of the sam¬ 
ple. We then individually consider the 112 objects (representing 
only 0.5% of the 21 065 selected galaxies) lying below the aver¬ 
age Dn(4000) versus W2-W3 relation by more than 3cr. In nine 
cases a bright star is located close to the galaxy’s center and is 
covered by the SDSS aperture; this causes the reduced Dn(4000) 
value. In 11 objects we are dealing with clear misidentifications, 
since they are all high-redshift QSOs (z = 1.5-2.3). 

There are still 92 objects. When looking at their SDSS spec¬ 
tra, it emerges that most of them have a pronounced absorption 
corresponding to the Balm er lines. These ob jects are often re¬ 
ferred to as Eh-A galaxies (iGoto et al.ll2003h . The presence of 
the A stars produces a substantial contribution in the 4000 A 
spectral region, leading to low Dn(4000) values. The lack of the 
corresponding increase in the WISE color is due to the spectral 
shape of the A stars, with a sharp drop below the Balmer break. 
The paucity of UV photons, which are most effectively absorbed 
by dust and are consequently the dominant mechanism of dust 
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Fig. 3. Dn(4000) versus W2-W3 diagram. 
Black filled circles represent the relationship 
derived for the MPA-JHU sample that is shown 
in Fig. [T] The blue line is a polynomial fit to 
the l-cr limit of this relationship. Gray dots rep¬ 
resent the sources in the BH12-A sample with 
z < 0.1 and S/N> 2 in the W2 and W3 bands. 
Red dots highlight the objects with rest-frame 
emission line EW< 5 A and a FIRST central 
counterpart. Red dots circled in black are the 
nine BL Lacs candidates, i.e. those sources with 
Dn(4000) < 1.735 (shown by the horizontal 
dashed line) that meet the additional constraint 
on the H6 index (see Fig. 0. The boundaries of 
the “forbidden” zone and ‘wedge” region are 
highlighted with a green and a cyan line, re¬ 
spectively. 


heating, accounts for the reduced MIR emission. Nonetheless, 
such sources can be easily recog nized by measuring their Hi5 in¬ 
dex (IWorthev & Ottavianilll99'^ . 

In Fig. |2] (left panel) we show the relation between the 
Dn(4000) and indices for the same sample of SDSS galax¬ 
ies as used above. The Hd index monotonously decreases with 
increasing Dn(4000). 

In Fig. |2] (right panel), we plot the 92 objects located in the 
“forbidden zone” in Fig. [T] in the same diagram. They clearly 
separate into two sequences. The first follow the general H<5 
versus Dn(4000) trend, and they are most likely E-i-A galax¬ 
ies. Their position within the sequence depends on the fractional 
contribution of A stars in each given object. The second group 
of objects is instead located where BL Lacs are expected. This 
region connects the bottom right hand corner of the diagram, 
which is characteristic of quiescent galaxies, to the location of a 
featureless source, having Dn(4000) = 1 and H5=0. Indeed, five 
out of ten of them are associated with a radio source with a flux 
density higher than ~ 50 mJy, which is a clear indication for the 
presence of non-thermal emission. We recover them as BL Lac 
candidates in Sect. |3] 

Summarizing, the proposed tool for isolating BL Lacs from a 
large population of galaxies based on the Dn(4000) versus WISE 
W2-W3 color diagram appears to be promising, particularly for 
HBL objects, which are expected to populate a “forbidden zone” 
in this diagram. However, the presence of outliers in this region 
requires 

i) visually exploring the SDSS data to remove objects whose 
spectrum is contaminated by a nearby star, 

ii) individually checking for misidentified objects, 

iii) most importantly, including a further diagram formed by the 
Dn(4000) and H<5 indices in the analysis in order to exclude 
L-hA galaxies. 

After such an analysis we are left with only five possible con¬ 
taminants to the BL Lacs population, i.e. five galaxies (out of 
21 065) with a low value of Dn(4000) but lack both MIR ex¬ 
cess and radio emission. The nature of these objects is unclear; 
however, due to their extreme rarity, we do not expect them to 
significantly affect the results of our study. 



Dn(4000) 

Fig. 4. index versus Dn(4000) diagram. Black filled circles 
represent the relationship derived for the MPA-JHU sample that 
is shown in Lig. |2] The blue line is a polynomial fit to the l-cr 
limit of this relationship. All objects above this line are con¬ 
sidered E-hA contaminants. Gray dots refer to the sources in 
the BH12-A sample with z < 0.1 and S/N> 2 in the W2 and 
W3 bands. Red dots highlight the objects with the emission line 
EW< 5 A and a LIRST central counterpart. Red dots circled in 
black are the nine BL Lacs candidates. The vertical dashed line 
indicates the limit Dn^ax = 1.735. 


3. Looking for the least luminous BL Lacs 

We looked for the least luminous BL Lacs by considering the 
sample of 18 286 radio sources built bv iBest & Heckm^(l2012l 
hereafter BH12). They combined the MPA-JHU sample with 
the NVSS and the LIRST surveys, the cross-matching going 
down to a flux density level of 5 mJy in the NVSS. We first 
focused on the 7302 objects they classified as radio-AGNs, sep¬ 
arating them from the star-forming galaxies. These authors also 
restricted their work to radio sources within the “main galaxy 
sample” (IStrauss et al.l l2003) . comprising those galaxies with 
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magnitudes in the range 14.5 < r < 17.77 and redshift range 
0.01 < z < 0.3. The 7302 sources (hereafter the BH12-A sam¬ 
ple) were cross-matched with the AllWISE catalog available at 
http://irsa.ipac.caltech.edu/ using a 3 arcsec search 
radius and finding 7242 counterpart^ 

Since our aim is to detect the faintest BL Lacs, we started 
to analyze the closest Universe, up to z = 0.1. In Lig. |3]we 
present the Dn(4000) versus W2-W3 plot for all objects with 
z < 0.1 and with S/N of WISE data in the W2 and W3 bands 
> 2. A comparison with Lig. [T] shows that most radio galaxies 
are located in the same area as the MPA-JHU sample and that 
we are now dealing with mostly weakly active galaxies, with 
Dn(4000) ~ 2 and W2-W3 in the 0.5-1.5 range, while the locus 
of star-forming galaxies is poorly populated. When further se¬ 
lecting objects with i) rest-frame EW of the main emission lines 
< 5 A and ii) existence of a LIRST central counterpart to exclude 
double sources lacking a radio core, we are left with 609 sources 
(Fig.0, only a few of which lying in the “forbidden zone”. 

We searched for BL Lac candidates among those selected 
objects whose ratio between the jet to galactic emission at 3900 
Ais/ > 1/3, corresponding to Dn = (Dn°-!-/)/(1 + f) < 1.735, 
having adopted a Dn(4000) index for quiescent galaxies Dn° = 
1.98 (see the previous sectionfl The choice of the / lower limit 
is arbitrary, since it is a compromise between the possibility of 
detecting a non-thermal signature (hence distinguishing a BL 
Lac candidate from a passive galaxy) and the requirement of 
completeness. However, this assumption will be included self- 
consistently in the analysis presented in the following sections. 
In Lig. [3]there are 12 objects with Dn(4000) < Dnmax = 1.735 
and at the same time below the 1 -cr limit traced by the MPA-JHU 
sample (see Sect. 0. These are labeled with an ID number. 

Ligure |4] shows the H^ versus Dn(4000) diagram. If we re¬ 
quire that H5 is lower than the 1-cr limit traced by the MPA-JHU 
sample to exclude contaminating E-hA galaxies (post-starburst) 
as explained in Sect.|2] we are left with nine sources. Their prop¬ 
erties are listed in Table[T] Six of them are knowr^BL Lacs with 
Dn(4000) ~ 1.3-1.4, while the other three are new candidates 
with Dn(4000) values of about 1.7. 

As an example in Lig. |5]we show a SED and SDSS spec¬ 
trum of one of the selected objects, namely the BL Lac candi¬ 
date ID=2440 at z = 0.076. The departure of the SDSS spectrum 
and WISE data from the properly normalized galactic templates 
highlights the reduced Dn(4000) value of the source, as well as 
the W3 (and W4) excess, we interpret as due to the presence of 
the jet emission. 

3.1. Expanding the search fo 0.1 < z < 0.15 

The BL Lac candidates we found in the z < 0.1 range have all 
radio flux densities greater than 18 mJy, which is well above 
the limit of the BH12 sample, 5 mJy. This means that we can 
detect similar objects at slightly higher redshifts. As a result, 
to improve the statistical significance of our work, we expand 
our search for BL Lac candidates to a three times larger space 
volume, i.e. to z = 0.15. We follow the same method as presented 
in the previous section, but we recalibrate all relationships with 
objects in the proper redshift range (see Lig. |6]l. 

^ Actually 7271, but 29 are double identifications of which we kept 
only the closest object 

Dn(4000)= 1.735 corresponds to C = 0.423. 

^ They are in cluded in the _ catalog of blazars 

“Roma-BZCAT” jMassaro et all I20091) . available at 
http://WWW.asdc.asi.it/bzcat/ 



log V [Hz] 



Fig. 5. Top: SED of a new BL Lac candidate, ID=2440. Red cir¬ 
cles display WISE and dereddened 2MASS flux densities. The 
SDSS dereddened spectrum is shown in blue. Gray lines rep¬ 
resent two SWIRE templates of a 13 Gyr old elliptical galaxy 
normalized to the W2 band (solid) or to the SDSS spectrum 
(dotted). The departure of the SDSS spectrum from the latter 
template highlights the reduced Dn(4000) value of the source, 
while the W3 (and W4) excess appears from the comparison of 
the WISE data with the first template. Bottom; rest-frame SDSS 
spectrum with indication of the main AGN emission lines (upper 
labels) and galactic absorption lines (lower labels). 


In the redshift range 0.1 < z < 0.15, the Dn(4000) limit that 
corresponds to / = 1/3 is 1.713, since the average Dn(4000) 
value of quiescent galaxies in this redshift bin is reduced to 1.95. 
Our selection criteria identify 742 objects with weak emission 
lines. The Dn(4000) versus W2-W3 diagram (Lig. |7| is very 
similar to the one obtained at lower redshift, but there are more 
sources with weak lines in the “forbidden zone”. In particular, 
there are ten objects with low Dn(4000), below 1.4, and well 
separated from the other sources. The same objects stand alone 
in the H5 versus Dn(4000) plot shown in Lig. [8] 

Moreover, there are 16 other objects that are below the 1-cr 
limit of the mean Dn(4000) versus W2-W3 relationship. Twelve 
of them are also below the l-cr limit of the mean H5 versus 
Dn(4000) relationship. The SDSS image of the source with ID 
3959 shows a contaminating star, so we discard it. The remain¬ 
ing 11 sources are BL Lac candidates. 

The properties of all candidates are reported in Table[T] 

3.2. Completing the analysis of the BH12 sample 

So far we only considered the 7302 objects included in the 
BH12-A radio-AGN sample, i.e. sources from the main galaxy 
sample where we excluded star-forming galaxies. We checked 
whether some BL Lacs could hide among the remaining BH12 
10 984 objects (hereafter BH12-B sample), 10747 of which have 
a WISE counterpart. If we consider all the sources with z < 0.15 
that meet the selection criteria discussed above, we are left with 
377 objects. These reduce to 13 if we require that the position 
in the Dn(4000) versus W2-W3 is below the average relation 
defined for the MPA-JHU sample. However, one source is con¬ 
taminated by a close star, two objects are QSO with mismatched 
redshift, one has an extremely uncertain redshift. Ligure|9]shows 
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Table 1. BL Lac objects and candidates in the BH12-A sample at z < 0.1 


ID 

R.A. 

Dec. 


Fi^gHz log Lradio 

BZCAT 

BH12-A sample with z < 0.1 

507 

118.654457 

39.17994 

0.096 

48.80 

40.20 

BZBJ0754+3910 

2061 

119.695808 

27.08766 

0.099 

63.50 

40.34 

BZBJ0758+2705 

1906 

122.412010 

34.92701 

0.082 

154.9 

40.56 

BZBJ0809+3455 

5997 

127.270111 

17.90440 

0.089 

231.8 

40.81 

BZBJ0829+1754 

786 

193.445877 

3.44177 

0.066 

78.50 

40.06 

BZBJ1253+0326 

4156 

196.580185 

11.22771 

0.086 

117.9 

40.48 

- 

4601 

199.005905 

8.58712 

0.051 

18.1 

39.19 

- 

2440 

212.955994 

52.81670 

0.076 

321.6 

40.81 

- 

4022 

233.009293 

30.27470 

0.065 

61.20 

39.95 

BZBJ1532+3016 

BH12-A sample with 0.1 < z < 0.15 

1315 

14.083661 

-9.60826 

0.103 

99.3 

40.58 

BZBJ0056-0936 

1805 

119.060333 

27.40993 

0.140 

20.2 

40.17 

- 

1903 

121.805740 

34.49784 

0.139 

32.2 

40.36 

- 

2083 

132.650848 

34.92296 

0.145 

31.2 

40.39 

BZBJ0850+3455 

2281 

144.188004 

5.15749 

0.131 

30.3 

40.29 

BZUJ0936+0509 

2341 

151.793533 

50.39902 

0.133 

30.6 

40.30 

BZBJ1007+5023 

6152 

162.411652 

27.70362 

0.144 

15.6 

40.08 

- 

1845 

163.433868 

49.49889 

0.140 

56.2 

40.62 

BZBJ1053+4929 

6428 

169.276047 

20.23538 

0.138 

117.5 

40.92 

BZBJl 117+2014 

3591 

172.926163 

47.00241 

0.126 

127.4 

40.87 

- 

3951 

183.795746 

7.53463 

0.136 

81.7 

40.75 

BZBJ1215+0732 

3958 

185.383606 

8.36228 

0.132 

150.1 

40.98 

BZBJ 1221+0821 

6943 

194.383087 

24.21118 

0.140 

10.6 

39.89 

BZBJ1257+2412 

1768 

201.445267 

5.41502 

0.135 

16.6 

40.05 

- 

7186 

211.548355 

22.31630 

0.128 

8.3 

39.70 

- 

6982 

212.616913 

14.64450 

0.144 

434.4 

41.53 

- 

3091 

216.876160 

54.15659 

0.106 

28.5 

40.06 

BZBJ1427+5409 

3403 

217.135864 

42.67252 

0.129 

43.6 

40.43 

BZBJ1428+4240 

7223 

223.784271 

19.33760 

0.115 

11.2 

39.73 

- 

4756 

229.690521 

6.23225 

0.102 

210.9 

40.89 

- 

1087 

234.945602 

3.47208 

0.131 

11.3 

39.86 

- 

BH12-B sample with z < 0.15 

1103 

164.657227 

56.469769 

0.143 

221.4 

41.23 

BZBJ1058+5628 

1936 

176.270935 

19.606350 

0.022 

545.4 

39.91 

- 

5076 

180.764603 

60.521980 

0.065 

171.4 

40.39 

BZBJ1203+6031 

9591 

184.054611 

49.380711 

0.146 

9.5 

39.88 

- 

9640 

227.671326 

33.584648 

0.114 

4.5 

39.33 

BZBJ1510+3335 

10537 

233.696716 

37.265148 

0.143 

22.1 

40.23 

BZBJ1534+3715 


Column description: 1) ID number, 2) R.A. [deg], 3) Dec. [deg], 4) redshift, 5) flux density at 1.4 GHz of the FIRST central component [mjy], 6) 
the corresponding luminosity [ erg s“*] having adopted a null spectral index for the k-correction, 7) BZCAT name. 
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Fig. 6. Same as Fig.[T](left panel) and Fig.|2(left panel), but for 0.1 < z < 0.15. 
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Fig. 7. Same as Fig.[3l but for 0.1 < z < 0.15. 


the Dn(4000) versus W2-W3 diagram for the BH12-B sample. 
The remaining nine sources are indicated. When considering that 
three of them (ID 2309, 8058, 8345) have high H(5 index and are 
thus likely E+A galaxy contaminants, we are left with six BL 
Lac candidates. They are listed in Table[T] 

We stress that the constraint based on the W2-W3 color al¬ 
lowed us to remove 54 objects with rest-frame EW of the main 
emission lines < 5A that would instead be selected as candi¬ 
dates based on the Dn(4000) index alone. Their location in the 
Dn(4000) versus W2-W3 diagram suggests that these objects 
are low-luminosity LINERs and/or Seyfert galaxies (see Fig.[TJ. 
Including the MIR data thus substantially improves the purity of 
the resulting BL Lac sample. 

3.3. Properties of the selected objects 

Our method allowed us to isolate 36 BL Lac objects (or can¬ 
didates). A detailed analysis of each individual source is de¬ 
ferred to a future paper. Nonetheless, we note that among the 
15 new BL Lac candidates, ten are associated with an X-ray 
source according to the ASDC SED builder too0, supporting 
their identification as genuine BL Lacs. Moreover, all but one 
object are massive elliptical galaxies, with stellar velocity dis- 
persion cr > 160km s ', which means togiM^n/Mo) > 7.8 
(iTremaine et al.l2()03) . i.e. the typical hosts of radio-loud AGNs. 
Their radio morphologies are consistent w ith those characteriz¬ 
ing BL Lacs (lAntonucci & Ulvestadlll985h : i.e., they are domi¬ 
nated by a compact radio core in the FIRST images, and in some 
cases, they also show a one-sided jet or a halo. 

We derive the radio-opticaQ spectral index of the non- 
thermal emission for our candidates by estimating the jet optical 
contribution from the observed Dn(4000) value; i.e. 

^ D“-Dn(4000)^ 

^i.o - 0 _ 1 

® http://tools.asdc.asi.it/SED/ 

’ Estimated between 1.4 GHz and 3900 A, both in the source rest- 
frame; we adopted a null radio spectral index and no k-correction is 
applied. 



Dn(4000) 


Fig. 8. Same as Fig. 0] but for 0.1 < z < 0.15. 


where Fo is the SDSS flux density (reddening-corrected) at 3900 
A. Figure [10] shows ctro versus the logarithm of the radio lumi¬ 
nosity of the FIRST central component L,. The mean value is 
< Pro >= 0.49 and the standard deviation 0.08. The range of a^o 
covered by our objects overlaps that of the BL Lac s analyzed by 
ISambranaet alJ (ll996h . |Padovani et alJ (l2003h . and iDonato et al.l 
(1200 Ih . going from ~ 0.2 to ~ 0.75. 

Figure [TT] shows the location of our BL Lac candidates in 
the WISE color-color diagram. They follow a sequence, start¬ 
ing f rom the region where passive elliptical galaxies are located 
(e.g.. lWright etal.ll20 lot [Raiteri et alj2014t) . and moving toward 
the so -called “WISE blazar strip” identified by iMassaro et aH 
(1201 ih . We also show tracks corresponding to the addition of a 
power-law emission component Fy ~ of increasing strength 
to a z = 0.1 elliptical galaxy flux. These tracks were obtained 
starting from the 13 Gyr elliptical galax y template of the SWIRE 
Template Librarj0 (iPolletta et al.ll200^ for different values of a 
in the range 0-1. They cover the region where most of our can- 


http://www.iasf-milano.inaf.it/~polletta/templates/swire_templates.html 
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Fig. 9. Same as Fig. [3] but for the BH12-B sam¬ 
ple at redshift z < 0.15. 


0.8 
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0.6 

0.5 
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0.3 

o.e 

39.0 39.5 40.0 40.5 41.0 41.5 

log 

Fig. 10. Radio-optical spectral index am versus the logarithm 
of the radio luminosity Ly for the 36 BL Lac candidates in the 
“forbidden” zone of the Dn(4000) versus W2-W3 diagram (blue 
dots). The red squares correspond to the five objects located in 
the “wedge” region of the same diagram. 



didates are found. This supports the presence of a power-law 
emission component in the MIR bands, whose relative weight 
increases at longer wavelengths. 

Overall, the properties of our candidates are reassuringly 
consistent with those of classical BL Lacs. 


4. Completeness of the low-luminosity BL Lacs 
sample 

Figure [12] shows the radio luminosity distribution of the FIRST 
central component, L,, for the 36 BL Lacs identified with our 
method. They extend up to log L, ~ 41.5 erg s *. The number of 
BL Lac candidates rapidly increases toward lower L,, reaching a 
peak at log Lj ~ 40.5 erg s“^ and then falls. This drop can be due 
to a genuine paucity of low-luminosity objects or, alternatively. 



Fig. 11. WISE color-color diagram for the BH12 sample with 
S/N(W3) > 2 and z < 0.1 (gray dots). Red hlled circles represent 
our BL Lac candidates in the same redshift range. The black 
tracks indicate the sequence followed by objects where a power- 
law component Fy ~ v““ of increasing strength is added to the 
emission of an elliptical galaxy, represented by the black dot. 
This was derived from the SWIRE template of a 13 Gyr elliptical 
galaxy at z = 0.1. Different tracks correspond to different values 
of a in the range 0-1. 


to a decrease in the completeness of our sample. Before deriving 
the BL Lac radio luminosity function (RLF), we must consider 
the selection effects included in the various steps described in 
Section|3] 

4.1. Completeness of the NVSS-SDSS sample 

From the point of view of the radio completeness, BH12 in¬ 
cluded the sources with a NVSS flux density larger than 5 mJy. 
The vast majority of our selected BL Lac candidates have fluxes 
that are much larger than this limit, an indication that the flux 
threshold does not have a signihcant impact on our results. 
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20 


15 


S: 10 


5 


38 39 40 41 42 43 

log i. 

Fig. 12. Radio luminosity distribution at 1.4 GHz of the 36 BL 
Lac candidates. The dotted black histogram includes the five 
sources in the “wedge”, while the green one also considers the 
two BL Lacs tagged as “STAR” in the SDSS, see Sects. 1431 and 



Fig. 13. Distribution of r band magnitudes of the galaxies in 
the BH12-A sample with 0.05 < z < 0.15. The vertical dotted 
lines indicate the limits defining the SDSS main galaxies sample. 
The red histogram reports the SDSS completeness in percentage 
from iMontero-Dorta & Pradj (l2009l) . 


Nonetheless, we include the radio flux limit in the simulations 
presented in Sect. |5] 

As for the optical s electio n of the sample, according to 
iMontero-Dorta & Pradal (l2009l) the redshift completeness of the 
SDSS decreases with decreasing apparent magnitude, starting 
from ~ 90% at the SDSS spectroscopic limit of r = 17.77 and 
reaching ~50% at r = 11.75 (see Fig. [nil. Most of the incom¬ 
pleteness is due to the SDSS fiber cladding that prevents fibers 
on any given plate from being placed closer than 55" apart. This 
does not introd uce a bias into the sample, but leads to a random 
loss of ~ 10% (IZehavi et al.ll2002ll of the potential spectroscopic 
targets. For brighter (and more extended) objects, other effects 
become important, such as the superposition of bright saturated 
stars on the target. 

Although the radio-loud AGN (RL-AGN) hosts are very lu¬ 
minous galaxies, it is important to assess their completeness 
level from the point of view of the SDSS spectroscopy. In Fig. 
[13] we show the distribution of r band magnitudes of the galax¬ 
ies in the BH12-A sample with 0.05 < z < 0.15. The distribution 
peaks at r ~ 16.2, and it rapidly falls on both sides. The num¬ 
ber of objects is essentially reduced to zero before reaching the 
high and low magnitude limits imposed by their selection in the 
SDSS main galaxies sample (r=17.77 and 14.5, respectively). 
This indicates that the selected sample of galaxies is complete. 

Conversely, very nearby (and bright) galaxies might fall in 
the region of low completeness at low magnitudes. We then pre¬ 
fer to set a more stringent lower limit of z = 0.05 on the galaxies 
considered. From the point of view of the BL Lac candidates, 
we discard only one of them (ID 1936), while from the point of 
view of the volume sampled, this is reduced by just 4%. 

4.2. Check on the WISE constraints 

Could our requirements on the existence of a WISE counterpart 
and/or on the quality of WISE data affect the BL Lac candidates 
count? Among the 60 sources in the BH12 AGN sample without 
a WISE counterpart, there are no objects meeting our criteria to 
be considered BL Lac candidates; EW less than 5 A, existence 
of a FIRST counterpart, and Dn < Dn^ax- 


If we consider the 445 sources with S/N(W3)< 2 in the 
BH12-A sample, only three objects meet the selection criteria. 
However, the low quality of the overall multiband dataset for 
these sources suggests that we should not include them among 
the BL Lac candidates. In fact, not only are the WISE data very 
uncertain, but the optical indices also have large errors. We ver¬ 
ified that among the 237 objects in the BH12-B sample without 
a WISE counterpart or with low S/N(W3), there are no BL Lac 
candidates. 

4.3. Are there BL Lacs outside the forbidden zone? 

The selection method we propose is based on the idea that the 
BL Lacs are found in the “forbidden” zone, while other radio 
emitting galaxies are located along a well-defined locus in the 
Dn(4000) vs. W2-W3 plane. However, there is the possibility 
that some BL Lacs follow this same sequence. Indeed, the sim¬ 
ulated tracks presented in the right hand panel of Fig.[T]indicate 
that the objects with a spectral index in the MIR - optical bands 
> 1 overlap with this region. 

We then consider those sources that meet the same require¬ 
ments as the ones used to select BL Lacs objects (rest-frame 
EW of the emission lines < 5 A and existence of a EIRST cen¬ 
tral counterpart) and that are located in the “wedge” region of 
the Dn(4000) vs. W2-W3 plane. This is the region below the 
Dn(4000) limit but above the polynomial fit to the 1 - cr limit of 
the relationship derived for the MPA-JHU sample (see Eig. [3]l. 
We And 54 sources, 22 of which have 0.05 < z < 0.15. When 
discarding the Eh- A galaxies, their number is reduced to ten. Eive 
of them are spiral galaxies, not the typical BL Lac hosts, and 
they are most likely low-luminosity AGNs. The remaining five 
are elliptical galaxies: while one is a BZCAT object (probably a 
genuine BL Lac), the nature of the other four is highly uncertain. 
Their compact radio morphology does not provide us with any 
useful information. We list these objects in Table |2l They have 
Dn(4000) values of 1.6-1.7 and W2 - W3 colors in the range 
1.7-2.1; i.e., they are all concentrated in a small portion of the 
“wedge” region, close to its apex. Their radio luminosity ranges 
from logLr ~ 39.5 to 40.3 [ ergs^'], see Eig. [121 their a^o 
goes from 0.44 to 0.57 (see Eig.fTOl). 
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Table 2. Possible BL Lacs in the BH12 sample with 0.05 < z < 0.15 in the “wedge” region of the Dn versus W2 - W3 diagram. 


ID 

R.A. 

Dec. 

W2-W3 

Dn(4000) 

Z 

Fi.4GHz 

l0§ ^radio 

BZCAT 

935 

147.13342 

55.59314 

2.064 

1.610 

0.118 

28.6 

40.16 

BZUJ0948-t5535 

3062 

194.97218 

57.86380 

1.711 

1.707 

0.149 

9.9 

39.92 

- 

5525 

213.01450 

29.46714 

1.970 

1.619 

0.115 

37.0 

40.25 

- 

2468 

226.09259 

47.68667 

1.957 

1.721 

0.093 

9.5 

39.46 

- 

1278 

359.67532 

-9.90999 

2.146 

1.663 

0.105 

31.8 

40.10 

- 


Column description: 1) ID number, 2) R.A. [deg], 3) Dec. [deg], 4) W2-W3, 5) Dn(4000), 6) redshift, 7) flux density at 1.4 GHz of the FIRST 
central component [mjy], 8) the corresponding luminosity [ erg s“'], 9) BZCAT name. 


4.4. An external completeness test: the BZCAT 

We here assess the completeness of our sample against the com¬ 
pilation of BL Lac objects (or candidates) included in the multi¬ 
frequency catalog of blazars Roma-BZCAT. The edition 4.1.1 
(August 2012) contains 3149 blazars, 1221 of which are BL 
Lac objects or candidates (BZB), and 221 are blazars of uncer¬ 
tain type (BZU). Such objects have been selected with different 
methods and from observations in various bands. The compar¬ 
ison might unveil the presence of bias in our selection method. 
Details on the cross-match between the BZCAT sources and our 
sample are given in Appendix lAl 

In summary, in the 0.05 < z < 0.1 redshift range, the com¬ 
parison between the BH12 sample and BZCAT suggests that we 
did not miss any object because of our selection method. We 
only lost 3 objects out of 18 because of fiber collision, which 
does not, however, introduce any bias into the sample but only a 
lower effective sky coverage. 

Moving to higher redshift, we recover all BZB and lineless 
BZU with the exception of six objects: two are bright featureless 
sources that are not included in the BH12 sample because they 
are tagged as “STAR”. We consider them in the RLF discussion. 
One object was instead missed because it lies in the “wedge” 
region of our Dn(4000) versus W2-W3 diagnostic plane. As 
shown in Fig. [1] this is the region populated by sources with 
MIR-optical spectral index a ~ 1, so likely LBL objects. Finally, 
three were discarded by our method because of a high Dn(4000) 
value, between 1.95 and 2.01. However, two of them also have 
WISE colors that are typical of quiescent galaxies, W2-W3 = 
0.76 and 0.39, respectively. For one of them, as detailed in the 
Appendix, we have strong evidence that it might be an incorrect 
identification. 


5. BL Lac selection function 

As discussed in Sect. [3 the main requirement for including an 
object belonging to the NVSS-SDSS sample in our candidates 
list is a flux contribution from the non-thermal jet component, 
Lj o, of at least one-third of the galactic one at 3900 A (rest 
frame). 

To estimate how this criterion is reflected in the selection 
function at different values of radio luminosity, we ran a Monte 
Carlo simulation. In each bin, we extracted 100 000 objects 
with: 

i) a random distribution in redshift between 0.05 < z < 0.15, 
properly weighting for the appropriate volumes and assum¬ 
ing a constant source density with redshift, to derive the cor¬ 
responding radio flux; objects with radio flux density below 
the 5 mJy threshold of the BH12 sample are discarded; 

ii) a random spectral index a^o (see below for a detailed 
discussion on this parameter) that leads to an estimate, 
for each object, of the expected jet optical emission 



Log Lj. 

Fig. 15. Expected fraction of BL Lac candidates recovered with 
our selection method at various levels of radio luminosity ob¬ 
tained with Monte Carlo simulations. The various lines refer to 
different assumptions on the radio-optical spectral index. Red: 
cTro = 0.49, the mean obtained from the 35 BL Lac candidates 
with 0.05 < z < 0.15; green: cTro = 0.53, the mean of the candi¬ 
dates with logLi > 40.5 ergs”*; blue: = 0.59, the mean of 

the IJy BL Lacs. 


[ erg s * cm ^ A * ] as 

logi^^ = logF, [mJy] - 5.74 - 14.71; 

iii) a random value of the flux density at 3900 A (rest frame) 
with average and dispersion measured from the sample 
galaxies divided in redshift bins of 0.01. In Fig.[T4]we show 
the distribution of radio versus optical fluxes into two red¬ 
shift ranges. 


For each object we tested whether the ratio between jet and 
galactic optical emission is above the f - 113 threshold. The 
fraction of objects expected to be selected as BL Lac candidates 
as a function of radio luminosity can Anally be derived. 

The assumption on the distribution of the radio-optical spec¬ 
tral index is the key to this simulation. In Sect. l3.3l we estimated 
a mean value of < a^o >= 0.49 + 0.08. However, our results 
might be plagued by selection effects, since high a^o objects of 
low radio luminosity could escape detection in the optical band. 
To assess how robust our estimates are, we compared them to 
those available in the literature. The avera ge spectral index of 
the (radio-selected) 1 Jy sample of BL Lacs (IStickel e t al.ll 19911) 
is < am >= 0.59 with a dispersion of 0.10 (ISambruna et al.l 
Il996h . Adopting their average for the simulation, we ob¬ 
tain the selection function reproduced as the blue curve in Fig. 
M it reaches 80% at logLr ~ 40.5 [ ergs *]. This means that. 
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Fig. 14. Left panel; the flux density of the FIRST central component, Fr, versus dereddened flux density within the SDSS fiber at 
3900 A, Fo, for the sources with 0.05 < z < 0.1. The candidate BL Lacs are plotted with red circles. Right panel: same as in the left 
panel, but for 0.1 < z < 0.15. 


at least above this threshold, the BL Lacs sample we selected 
is very close to completeness and that the distribution of cHro 
is not significantly affected by selections biases. Its average, 
< am >= 0.53 + 0.08, can be considered as a robust value. The 
selection function obtained by adopting this value is shown in 

FigM 

The 1 Jy objects are, in general, much brighter than the BL 
Lacs we are considering, with radio luminosities reaching ~ 
10^^ erg s“'. If we limit to the 1 Jy sources with < lO'^^ erg s“' 
we obtain < am >- 0.52+0.09, which agrees remarkably with 
the value quoted above, while the average am for those with 
Lr > 10« ergs * is 0.65+0.08. These results suggest there is 
a trend between am and the rad io luminosity, as ex pected from 
the so-called “blazar sequence” (iFossati et al.ll99^ . This would 
account for the (slightly) lower average value we find for our 
LPBL sample. Indeed, the sources of the 1 Jy sample (limiting to 
the 22 with measured redshift) show a linear correlation between 
am and L^, with a slope of 0.06. The Spearman rank correlation 
test returns a probability of only 2% that the two variables are 
not correlated. We therefore repeat the simulation by adopting 
the observed distribution of spectral indices; i.e., < am >- 0.49 
(Fig.m. The selection function obtained with this assumption 
differs only marginally from what is derived for < am >- 0.53, 
with a horizontal shift of only ~0.2 dex. 

6. The BL Lacs luminosity function 

Several studies have dealt with the determination of the BL Lacs 
RLF In all cases, the approach was based on flux-limited sam¬ 
ples. In our study, instead, we consider an essentially volume- 
limited sample, because the radio, optical, and MIR flux lim¬ 
its do not represent a significant constraint in our selection. The 
advantage of this method is twofold; all objects have a known 
redshift, and they are all within a very limited range, so that evo¬ 
lution can be neglected. 

We can then derive the BL Lacs RLF by simply dividing 
the number of detected objects by the total volume covered, re¬ 
calling that the intersection between the spectroscopic DR7 and 
FIRST covers 17.3% of the whole sky. By also considering the 
effects of fiber cladding, this results in a volume of 0.158 Gpc^ 
from z = 0.05 to z = 0.15. We convert the radio luminosity 
distribution shown in Fig. [12] (including the 35 candidates with 


log L, [erg s“'] 


39.0 39.5 40.0 40.5 41.0 41.5 42.0 



log P, [W Hz"‘] 


Fig. 16. RLF at 5 GHz of the selected BL Lac candidates (black 
histogram). The gray histogram shows the effects of includ- 
ing the five ob j ects l ocated in the “wedge”. The results from 
iPadovani et al.l ( 2 0071) are in red, wh ile blue dots refer to the 
work by iMarcha & Caccianigal (1201 3h for “classical BL Lacs”. 
We also report (as plus sig ns) three points at low luminosity of 
the iMarchg & Caccianiga RLF for “Type 0” AGNs. The green 
line is the RLF normalized to match the number of objects of 
our sample with logPr > 23.8 [WHz *L afte r having adopted 
the same slope of iMarcha & Caccianigal (1201 3h . 


0.05 < z < 0.15 and the two objects tagged as “STAR”) by divid¬ 
ing the number of objects by the radio luminosity, the selection 
function, and the total volume. The errors on the RLF have been 
estimated by combining Poisson statistics with the uncertainties 
related to the selection functioifl The result is shown in Fig. [T6| 
To ease the comparison with previous studies, we present the ra- 

^ We considered the two assumptions: < > = 0.49 and 0.53. The 

differences in the resulting RLF are of only 0.17, 0.15, and 0.05 dex for 
the three bins at the lowest luminosities. By adopting instead < > 

= 0.59, these hins are further increased hy 0.75, 0.38, and 0.12 dex, 
respectively. 
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log L, [erg s '] 


39 40 41 42 



Fig. 17. Observed number of BL Lacs (black solid line) versus 
radio power. The blue histogram represents the predicted num¬ 
ber of objects in our volume extrapolating the RLF derived for 
our sample to low radio luminosities (green line in Fig. fTbl) and 
multiplied by the selection function (red line in Fig. fTsT l. The 
gray histogram shows the effects of the inclusion of the five 
objects located in the “wedge”. The red histogram displays the 
same quantity as for a RLF with B - 1.75 and with a break at 
log Pi - 24.5 [W Hz '], below which the RLF slope is B1 =0. 


dio luminosities in both cgs and MKS units and a ssume a zero 
radio spectral index to convert from 1.4 to 5 GHz dStickel et al.l 

[Mil). 

We can compare our results on the RLF with those ob¬ 
tained by pr evious works and, in p a rticula r, with iPadovani et al.l 
(l2007h and iMarcha & Caccianied (l2013h . We assume a RLF 
in the form of a power law, OjPr) oc Since our sample 
covers a limited range of radio luminosities, we refrain from 
calculating the R LF slope. If we adopt a pow er-law slope of 
B - 2.31 as in iMarcha & Caccianigal (1201 3h for “classical” 
BL Lacs, we can calculate the normalization of our RLF so 
that the integrated RLF returns the number of observed objects 
above logPr = 23.8[WHz '], the range of po wers common 
to both studies. This is displayed in Fig. llafT The compari- 
son shows very good agreement with the results obtained by 
IMarcha & Caccianigal for BL Lacs and a mismatch between our 
RLF and the Padovani et al.l one by a factor ~ 3. 

The main aim of our study is to isolate and analyze the 
least luminous BL Lacs. Indeed, we detected candidates down to 
log Pi ~ 23 [W Hz^']. We defer a detailed study of the individual 
objects and of their SED to a future paper. Here we want to ad¬ 
dress the behavior of their RLF in a poorly studied regime of low 
radio power. The only study r eaching these luminosities is from 
IMarcha & Caccianigal (l2013l) : they consider not only classical 
BL Lacs but also the objects they defined as “Type 0” and “weak- 
lined” AGNs. Their space density is substantially higher than 
for our BL Lac sample, with, say, about ten times more Type 0 
AGNs than BL Lacs between logP, = 23.5 and 24.0 [W Hz ']. 
However, recent studies have shown that the dominant popula¬ 
tion of radio galaxies in the local Universe is formed by the so- 


Had we used B = 2.12 as derived bv IPadovani et all ll2007tt . 
we would have obtained a difference of only 0.08 dex at logP^ = 
24.5 [WHz-1] 


called FR 0 sources (iBaldi & Capettil 120101 : ISadler et alJl20l4 
iBaldi et ^l20L5ll . Their lines EW are generally low (more than 
80% showing values lower than 5A, R. Baldi, private commu¬ 
nication). Moreover, they are generally sources with a high core 
dominance and, for this reason, ofte n flat spectrum radio sources. 
As explained by IBaldi et ^ (l20L5l) the deficit of extended radio 
emission is an intrinsic property of ER 0, and it is not due to a 
beaming effect of orientation, as in the case of BL Lacs. In other 
words, a signihcant fraction of the ER 0 population conforms to 
the definition of Type 0 and weak-lined AGNs. Therefore, the as¬ 
sociation of all Type 0 and weak-lined AGNs with LPBL might 
not be justified. 

We estimate the predicted number of LPBL in our volume 
by extrapolating the RLE derived above to low powers, taking 
our selection function into account (see Eig. dsi). The result is 
shown in Eig. [17] where we compare the observed and predicted 
number counts. A large disagreement emerges, where the differ¬ 
ence is a factor of ~ 10 already at log Pi ~ 10^^'^^ WHz ' (i.e., 
logLr ~ 40), where our selection function is high (around 70%) 
enough to make this result robust. 

We then consider the possibility of a break in the RLE occur¬ 
ring at Pbreak- Eurthermore, it appears that also at high luminosi¬ 
ties the slope might be less than the value of 2.31 adopted. We 
are not in the position of accurately determining all parameters 
describing the RLE with the available data. Nonetheless, in Eig. 
[m we show an illustrative example of a RLE with B - 1.75, a 
break at logPbieak = 24.5 [WHz“^], and a slope below the break 
51=0, which fairly matches the observed counts at all lumi¬ 
nosities. This result is effectively independent of the assumption 
used for the radio-optical spectral index; had we used the av¬ 
erage value derived for the 1 Jy sample, the break luminosity 
would have decreased by just a factor ~ 1.7. 

Values of 51 below 1, as suggested by this analysis, corre¬ 
spond to BL Lacs number counts growing with luminosity up to 
Pbreak, implying the existence of a peak. 


6.1. The BL Lac RLF and the AGN unified modei. 


According to the AGN unihed scheme, the RLEs of BL Lacs 
and RL-AGNs are expected to be linked to each other. In fact, 
blazars are th ought to be the beamed version of radio galax¬ 
ies (see, e.g.. lUrrv & Padovanilll99-^ . The effect of relativis¬ 
tic beaming on the RLE was analyzed in a series of papers 
by lUmv & Shafed ([I9PI). lUiry & Padovanil (Il991l) . lUrrv et all 
(Il99lh . and llJrrv & Padovanil (Il995 ). If the RLE of the parent 
population of unbeamed objects is described by a broken power 
law, the break will produce a change of slope in the BL Lac RLE 
at a power amplihed by ~ where is the maximum 

value of the jet Doppler factor distributioiO- 

The RLE of RL-AGNs extending to the lowest value of ra - 
dio luminosity has been obtained by iMauch & Sadlej (l2007h . 
covering the range logP^ = 20.4 - 26.4[WHz ']. At low lu¬ 
minosities, the RLE does not significantly depart from a power 
law, and their results suggest that any break must be located be¬ 
low logPr ~ 2L6WHz where the limited number of ob¬ 
served objects makes the RLE shape uncertain. On the other 
hand, the number density of radio emitting AGNs cannot ex- 
ceed that of their potenti al hosts: this argument has been used by 
ICanetti & Raiteril (l2015h . requiring a minimum black hole mass 


** The Doppler factor (5 is defined as 5 = [r(l - p cos0)] where 
r = (1 - is the hulk Lorentz factor, p is the velocity in units of 

the speed of light, and 6 the viewing angle. The exact value of the S 
exponent depends on the geometry of the emitting region. 
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of log (Mbh/ALq) > 7.5 for the host of RL-AGNs. With this ap¬ 
proach they found that the break in the RLF of RL-AGNs must 
be located at ^ 10^** ^ W Hz '. 

The ratio between the break powers in the RTFs of BL Lacs 
and RL-AGNs is a factor ~ 10^ - 10"^. If adopting p - 3, this 
suggests a value of the Doppler factor dmax ~ 10 - 20, consis¬ 
tent with the typi cal results from the observations and the mod el 
predictions (e.g.. lGhisellini et^ll993l:ISavolainen et al.ll2010l) . 


7. Summary and conclusions 

In this paper we presented a new strategy for recognizing BL 
Lac objects, based on the combination of optical indices de¬ 
rived from the SDSS spectra with MIR colors from the AllWISE 
survey. This enabled us to separate BL Lac objects from low- 
luminosity AGNs and star-forming galaxies. Our method is par¬ 
ticularly useful for identifying LPBL, where the traditional op¬ 
tical non-thermal signature is reduced thanks to the dominant 
host-galaxy emission. We applied it to the sample of 18 286 ra¬ 
dio sources selected by BH12 by combing the SDSS DR 7 spec¬ 
tra with the NVSS and FIRST surveys. We showed that in the 
Dn(4000) versus W2-W3 plane, BL Lacs fall in a region (de- 
hned as “forbidden zone”) that is scarcely populated by the gen¬ 
eral galaxies population. The main contaminants are the Ah-E 
(post-starburst) galaxies that can be easily removed with the 
aid of their characteristic pronounced Balmer absorption lines. 
When considering low luminosities, the use of the Dn(4000) 
index alone is not sufficient to obtain a clean sample. In fact, 
the constraint based on the W2-W3 color allowed us to remove 
54 objects (most of which are low-luminosity LINERs and/or 
Seyfert galaxies) that would instead be selected as BL Lac can¬ 
didates. Including the MIR data thus substantially improves the 
quality of the resulting BL Lac sample. 

We isolated 36 BL Lac candidates up to z = 0.15. All but 
one are associated with luminous elliptical galaxies with radio 
powers ranging from logLr = 39.2 to 41.5 [ ergs^']. Eifteen of 
the selected objects are new identihcations, most of which have 
an X-ray counterpart. 

We discussed various issues concerning our sample com¬ 
pleteness. The radio flux threshold does not have a signihcant 
impact on our results. As for the optical selection, most of the 
expected RL-AGN hosts are included in the MPA-JHU catalog. 
The only exceptions are the low-luminosity tail close to the up¬ 
per redshift boundary and the most luminous galaxies located at 
very low redshift. We then set a more stringent lower limit of 
z = 0.05 on the galaxies considered. We used the BZCAT for 
an external check of our results. We recovered all the known ob¬ 
jects in the catalog up to z = 0.1, while we lost less than 10% 
of them in the 0.1 < z < 0.15 range. In particular, only two 
very bright BL Lacs are lost because of their point-like appear¬ 
ance, which caused these objects to be misclassihed as “STAR” 
in the SDSS. We also tested the effects of the criterion based on 
the W2-W3 threshold. Out of the 54 sources mentioned above, 
we hnd five objects that conform to the general properties of 
BL Lacs (since hosted in elliptical galaxies and having low lines 
EW). However, they are of uncertain nature because they lie in 
a region of the W2-W3 vs. Dn(4000) plane populated by other 
classes of AGNs. 

We simulated the selection function of BL Lacs within our 
sample. The key ingredient is the assumption on the distribution 
of the radio-to-optical spectral index for which we adopted the 
values derived from our candidates. Since they agree with those 
found in the literature, our results are not signihcantly affected 


by selections biases. The sample completeness reaches ~ 50% 
at logLr ~ 39.7 [ergs^']. 

The selected BL Lac candidates were used to explore the BL 
Lacs RLF at low radio powers. Since we are dealing with an es¬ 
sentially volume-limited sample at low redshift, we can neglect 
evolutionary effects on the RLF. This can be obtained simply by 
counting all objects at the various radio luminosities within the 
considered volume. We found good agreement with the results 
by previous authors for log > 24.5 [W Hz ']. 

However, we found a dramatic paucity of LPBL with respect 
to the extrapolation of the BL Lacs RLF, a result that is altered 
neither by the inclusion of the hve potential BL Lacs nor by the 
uncertainties on the selection function. This occurs well above 
the power at which our selection function falls signihcantly and 
requires a break in the RLF. We estimated that the break is lo¬ 
cated at log Pr ~ 24.5 [W Hz ']. The slope of the RLF below the 
break is consistent with a zero value, implying the presence of a 
peak in the number density distribution of BL Lacs. 

In the framework of the AGN unihcation scheme, BL Lacs 
are the beamed counterparts of low-power RL-AGNs, so that a 
link between their RLFs is expected. In particular, a low-power 
break in the parent population translates into a change in slope 
in the RLF of BL Lacs. Indeed, there is evidence of a low- 
luminosity break in the RL-AGN RLF located about three to 
four orders of magnitude below what we found for BL Lacs. 
This corresponds to a maximum Doppler factor of ~ 10 - 20, in 
agreement with the observations and model predictions. 

In the near future, we plan to do a detailed analysis of the 
whole BL Lac candidate sample derived in this paper with two 
aims: i) to conhrm their nature as genuine BL Lacs and ii) to 
characterize their broad-band properties. In particular, we need 
to obtain high-resolution radio data in order to isolate their nu¬ 
clear emission and X-ray data to explore the high-energy non- 
thermal emission. 
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Nine of them are in our survey area, but two lack spectra because 
of fiber collision. The other seven have strong emission lines. 

In the redshift range 0.1 < z < 0.15, we found 59 BZB in 
BZCAT, 26 of which are located in the area covered by both the 
spectroscopic DR7 and FIRST. When excluding three objects 
that have no spectrum because of fiber collision and one object 
whose spectrum has a defect in the Ha region, we are left with 
22 sources, 15 of which belong to the BH12 sample. The other 
seven are lost because: 

- one object has no spectrum because it is too faint, but we know 
from BOSS that it is a spiral galaxy; 

- three are featureless objects whose SDSS redshift is higher than 
0.15; 

- one is a strong-lined object reaching an EW of ~ 10 A; 

- two (BZBJ0809-H5218 and BZBJ1221-H2813-W Comae) are 
classified as “STARS” in the SDSS and thus do not enter in the 
main galaxy sample analyzed by BH12; they have featureless 
spectra, and log =41.11 and 41.54, respectively. 

Among the 15 sources included in the BH12 sample, 13 
are in our list of BL Lacs candidates, while two sources 
(BZBJ081 Oh- 49 11 and BZBJ1404 h- 2701) were not selected by 
our method because of a high Dn(4000). They have logLr = 
39.50 and 40.39, respectively. The first one was identified as a 
BL Lac because of y-ray emission detected by the EGRET in¬ 
strument on-board the Compton Gamma-Ray Observatory, but 
it is not present in the third Lermi Large A rea Telescope source 
catalog. Lurthermore. iLiuzzo et al.l (1201 3h did not detect it with 
the VLA at 8.4 GHz, so we consider it a dubious association. 

The above analysis leads to the conclusion that in the 0.1 < 
z < 0.15 redshift range, our selection method recovered all BZB 
but two sources characterized by a high Dn(4000) value. Another 
two objects were lost because they are so bright that they are 
not classified as galaxies and so do not enter the BH12 sample. 
Another four objects lack complete spectral information. 

Among the BZU, there are 23 objects in BZCAT at 0.1 < z < 
0.15, 14 that are in the region covered by both the spectroscopic 
DR7 and LIRST. One object has no spectrum because of fiber 
collision, two have strong emission lines, 11 are in the BH12 
sample, but only three have rest-frame EWs less than 5 A and 
only one, BZUJ0936-I-0509, is selected by our method (Table 
□. The other two are BZUJ0820H-4853 and BZUJ0948H-5535. 
The former was discarded because of a high Dn(4000) of 2.01, 
while the latter lies in the “wedge” region of Lig. 7. They have 
logLr = 40.65 and 40.12, respectively. 

The analysis of the BZU in the 0.1<z<0.15 redshift range 
thus reveals that our methods did not recover two sources out of 
14. One other object was missed because of fiber collision. 


Appendix A: Details of the comparison with BZCAT. 

In the BZCAT list there are 34 BZB with 0.05 < z < 0.1, nine 
of which are located in the area covered by both LIRST and the 
spectroscopic DR 7. Seven of them are in our list of BL Lacs 
candidates. The remaining two are: 

- BZBJ1216-I-0929, which was selected as spectroscopic target, 
but its spectrum was not obtained owing to a fiber collision; and 

- BZBJ1426-1-2415, which shows a featureless spectrum; the lack 
of any signature from the galaxy emission cast doubts on the 
accuracy of its tentative redshift estimate, z = 0.055; indeed, the 
DR 7 value is z = 2.243. 

A check on the BZU is more problematic, because most of 
them are LSRQs. There are 24 BZU in BZCAT at 0.05 < z < 0.1. 
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